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We revisit the question of how a definite phase between Bose-Einstein condensates can sponta- 
neously appear under the effect of measurements. 

We first consider a system that is the juxtaposition of two subsystems in Fock states with high 
populations, and assume that successive individual position measurements are performed. Initially, 
the relative phase is totally undefined, and no interference effect takes place in the first position 
measurement. But, while successive measurements are accumulated, the relative phase becomes 
better and better defined, and a clear interference pattern emerges. It turns out that all observed 
results can be interpreted in terms of a pre-existing, but totally unknown, relative phase, which 
remains exactly constant during the experiment. 

We then generalize the results to more condensates. We also consider other initial quantum states 
than pure Fock states, and distinguish between intrinsic phase of a quantum state and phase in- 
duced by measurements. Finally, we examine the case of multiple condensates of spin states. We 
discuss a curious quantum effect, where the measurement of the spin angular momentum of a small 
number of particles can induce a big angular momentum in a much larger assembly of particles, 
even at an arbitrary distance. This spin observable can be macroscopic, analogous to the pointer of 
a measurement apparatus, which illustrates the non-locality of standard quantum mechanics with 
particular clarity. The effect can be described as the teleportation at arbitrary distances of the 
continuous classical result of a local experiment. The EPR argument, transposed to this case, takes 
a particularly convincing form since it does not involve incompatible measurements and deals only 
with macroscopic variables. 

PACS 03.65. Ta and Ud ; 03.75.Gg 

It is well known in quantum mechanics that one cannot construct a quantum state where the number of particles 
and the phase are both arbitrarily well defined; they are actually incompatible observables, related by a Heisenberg 
type uncertainty relation, as position and momentum of a single particle. For a general discussion of the phase 
operator in quantum mechanics, see for instance refs. An usual illustration of the relation between phase and 

particle number is given by the so called "coherent states" or "Schrodingcr semi-classical states" , where fluctuations 
of the number of particles are used in order to define a phase. These states are often discussed in the context of 
electromagnetism |4|, but they also apply to Bose-Einstein condensates, on which we focus the interest here - more 
precisely, multiple condensates and their relative phase. 

For instance we consider a physical system that is the juxtaposition of two systems in Fock states (number states), 
described by the state vector: 

\<Po> = iNa-.ipa ; Nb:ipb> (1) 

where Na particles are condensed in the same single particle quantum state ipa , and Nh particles in the same quantum 
state tpb- These two states may correspond to plane waves with momenta and k^, a case in which one can expect 
the occurrence of an interference pattern with spatial frequency k^ — k;,. Nevertheless, the position of this pattern 
depends of the relative phase of the waves, which is completely unknown in a state such as (Q. One could then wonder 
whether an interference effect is observable at all under these conditions. 

The question was studied theoretically by several authors |Elli0B@; see also refs. The result 

is that an interference pattern is in fact observed, with a phase that spontaneously emerges from the quantum 
measurement process itself. The detection of the first particle occurs at a completely random position, but this 
position provides a first information on the value of the phase; the second particle is then detected at a position that 
is correlated with the first, and this measurements makes the information on the phase more precise; as more and 
more detection events are accumulated, the phase becomes better and better defined. In practice, a relatively small 
number of detections is sufficient to determine the relative phase with good accuracy. Nevertheless, before the first 
measurement, according to standard quantum mechanics there is no way whatsoever to predict the value that this 
phase will spontaneously choose. In other words, repeating the same experiment from the same initial conditions 
will lead to another, completely independent, value of the relative phase. The spontaneous appearance of a relative 
phase under the effect of measurement is related to the notion of spontaneous symmetry breaking of the conservation 
of the number of particles, and to non-zero average values of the field operator in superfluid systems It also 

relates to Anderson's intriguing question, as quoted by Leggett ■ "Do superfluids which have never seen each other 
possess a definite relative phase?" [13, 13 • 
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In this article, we revisit this question by using a shghtly different method from that of refs. |H Q, |M |^ Hol ITll 
[T^ : we propose a straightforward analytical calculation that does not require the introduction of any (incoherent) 
fluctuation of the particle numbers, as in Q and for instance, but deals directly with pure Fock states; it reproduces 
the stochastic simulations of ^ with only one approximation (total number of particles much larger than the number 
of measurements). In f^jwe discuss how the process of measurement, applied to a pair of systems that are initially 
in Fock states, therefore with no well defined phase, can make such a phase appear spontaneously ; in f^we apply 
our calculations to more general cases and discuss how an intrinsic phase, contained in the initial state, can combine 
with a spontaneous phase induced by measurements; in illlll we discuss multiple condensates in the context of spin 
states and show in particular how the measurement of a small number of spins may induce an angular momentum for 
a much larger assembly of spins, a curious non-local quantum effect. 

I. SPONTANEOUS PHASE INDUCED BY MEASUREMENTS ON FOCK STATES 

We first recall a general result of quantum mechanics. We consider a system in a quantum state | ^'q >, and a series 
of observables: 

A with eigenvalues ai; P{ai) is the projector onto the corresponding eigenstate(s) 

B with eigenvalues hj] P{bj) is the projector onto the corresponding eigenstate(s) (2) 

C etc. etc. 

If we assume that all these observables commute with each other, then the probability of observing in a combined 
measurement the results at, bj, Ck, etc. can simply be written: 

< $0 I Pia,) P{bj) P{ck).. I $0 > (3) 

(the usual intrication of projectors in the so called Wigner formula is not necessary here because the operators 
commute). Whether or not the wave packet reduction postulate is applied does no matter; in other words, all 
measurements can be made simultaneously, or one after the other, and the order of measurements is irrelevant. What 
is assumed, nevertheless, is that the sequence of measurements covers a time that is negligible in comparison with the 
time constants associated with the intrinsic evolution of the physical system; this is why no evolution operator has to 
be introduced in the formula. We now apply lO to position measurements inside the overlap regions of two or more 
condensates. 



A. A simple case: two highly populated states 

We begin with a simple case, where the system is just a juxtaposition of two Bose-Einstein condensates in plane 
waves: 

I $0 > = I iVa : ; iVb : kb > (4) 

Here ka and k^ are single particle states of well defined momentum normalized in a box with periodic boundary 
conditions; we assume that Na and Nb are large numbers. The probability for detecting a particle at point r corresponds 
to the following operator, similar to the projectors P{ai) introduced above: 

N 

^^(r) ^'(r) = ^ I I : r >< i : r I (5) 

i=l 

Here ^'(r) is the usual field operator, and | i : r > is the one particle state corresponding to a perfect localization of 
particle i at point r; N = Na + Nt is the total number of particles. To make the operator \l/^'(r) Vl/(r) really similar to a 
projector with eigenvalues and 1, one has to integrate it over some small r domain Ar centered around r; assuming 
that Ar is sufficiently small ensures that the probability of finding two particles or more in Ar is negligible. For the 
sake of simplicity we do not write these integrations explicitly; in other words we write probability densities instead 
of probabilities, but it would be easy to come back to real probabilities by multiplying by an appropriate power of 
Ar. The field operator can be expanded onto the annihilation operators at of momentum states according to: 



^'(r) - ^'^"^ 

k 



(6) 
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where, as usual, the sum over k ranges over the values allowed by periodic boundary conditions in a box. 

We now calculate the probability for observing one particle at position ri, another at position r2, still another at 
position r3, .. one at position rp. We consider that all positions ri, r2, ..are different, so that all operators \[' and 4"^ 
commute (more precisely, we assume that the small integration domains An, Apji, .• do not overlap). This allows us 
to write this probability as: 

< $0 I *^(ri)*t(r2)«'t(r3)...^(r3)*(r2)*(ri) | % > (7) 

or, if is used: 

< $0 I a' ' a' ' a], a^p a^^ai^^ \ $o > 

Acting on I <i>o > given by the annihilation operators a's in the second line of this formula give zero unless the 
values of the k and k momenta are all equal to ka or k^. In addition, the series of k's must contain the same number 
of ka's (and k^'s) as the series of k 's; otherwise, the product of a and af operators creates a ket orthogonal to | $0 >, 
and the second line vanishes. When this condition is met, each Ok^ ^ or aj^ ^ introduces a factor ^ Na.b ± n, where n 
is some integer ranging from and P (the number of measurements). 

At this point we make our only approximation: we assume that P is much smaller than Na and iVf,, so that we can 
replace all these factors by yjNa.b- Now, in the first line of 0, consider for instance the two first summations over 
ki and kj^; if one of these vectors is chosen as k^,, the other as k^, one obtains a contribution ^/ NaN^e^^^^'^^^''^'^^ ; if 
both are equal to the same value, k^ or k^, one gets a contribution Na + Nf, = N . The same result is obviously also 
valid for all other sums. Therefore, if we introduce the notation: 

Fo(r) = N a + Nb = N 
where x is the contrast ratio x < x <1) defined by: 



we can rewrite ^ as: 



^ F,,(ri)F,,(r2) F,,{vp) (11) 



In this expression, the sum contains all possible values 0, ±1 of gi, ^2, etc., with the only condition that their sum 
be zero. Now, since the integral: 

(where g is a positive or negative integer) vanishes except if q = 0, we can release this condition by rewriting Ijllf) as: 

- / ■Y\{{Fo{v^) + e'''F+,{v,)+e-''•F^^{v.,)] (13) 

Finally, inserting Q into this result provides the probability in the form: 

- 1^ n {1 + ^ cos [(k„ - kfa) • r, + $]} (14) 

This result is simple to interpret: the dependence of any sequence of events is exactly the same as for a system 
with a well defined, but completely unknown, initial phase <&. For each value of this phase, the dependence of the 
probability factorizes and the detection events are independent; it gives exactly the usual interference pattern between 
two waves with wave vectors k^ and kf,. But the summation of $ between and 2tt destroys the factorization and 
correlates the events. 
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Let us assume for the sake of simplicity that the contrast ratio x is equal to 1. When the first detection localizes a 
particle at point Ri, some information on $ is immediately obtained, if only because the value $ — (k^ — k;,) • Ri + tt 
becomes incompatible with the first detection. After this measurement, the probability of a second detection at point 
V2 is proportional to: 

/ — (1+C0s[(k, -kfc)-Ri+$]){l+C0s[(k,-kb)T2+$]} (15) 

where the distribution of the values of $ is no longer uniform: it is now given by the Ri dependent sinusoidal function 
between the first brackets inside the integral. If two measurements are initially performed with localizations at points 
Ri and R2, the $ distribution is then given by the product of two sinusoidal functions, which is more peaked than the 
previous distribution, etc.. More and more precise information on the value of $ is progressively accumulated while 
more and more measurements are performed. At some point, the phase is practically determined and the detection 
events become quasi independent. As pointed out by M0lmer [l9|, this is because the two condensates tend to fuse 
into a single condensate under the effect of successive interference measurements. See also ref. for discussion of 
the evolution of the system towards a coherent state, which requires a number of detections comparable to the total 
number of atoms. For more details about the evolution of the $ distribution function, see for instance [tLIsLITiII. 

It is interesting to note that the phase plays a role which appears, mathematically, very similar to that of the 
so called "additional variables" (or "hidden variables") sometimes introduced to interpret the results of quantum 
mechanics - see for instance and for a review. If standard quantum mechanics is completed with an additional 
phase variable, one can consider that the measurement process reveals a pre-existing value of instead of creating 
the relative phase between the two condensates as in the standard interpretation of quantum mechanics. In our 
calculation, the motivation for introducing $ was not a fundamental rc-intcrpretation of quantum mechanics; it was 
just a convenient way to sum the many terms in © while maintaining the sum rule over the series of values of k and 
k . The precise role of the $ integral is to ensure that the successive creation and annihilation operators bring back 
the system to the same initial Fock state, in other words to enforce a number conservation rule. Each term in the sum 
can be seen as arising from one possible path over intermediates states in which the values of the occupation numbers 
vary in a given way before coming bask to their initial values. It is interesting to remark that the interference of all 
these paths should result in a simple phase integral, leading naturally to the introduction of an additional variable 
<&. We recover the usual relation between conjugate variables, where a sum over one variable with equal weights 
ensures a strict conservation of the other; the phase appears here as related to fluctuations of the number of particles 
in intermediate states (instead of in the initial state as usual). 



B. Three highly populated states 



by: 



We now generalize the preceding discussion to three states by assuming that the initial state of the system is given 



$0 > = \Na: ka; Ni, : k^; iV, : k^ > 



(16) 

etc. by using Q again. 



and evaluate the probability of particle position measurements at different points ri, r2, 
The calculation remains very similar to that of the preceding section. The only difference is that the k and k can 
now take three values, k^, kf, and kc, so that in the second line of (jSJl it is no longer sufficient to ensure that Na 
returns to the same value; the value of Ni, (or Nc) must also be controlled. This can be obtained by introducing two 
relative phases, $ and $ . Another difference is that the sums in the first line of JHJ now include three terms that 
contain the same k's, as well as six others that contain different k's and introduce cosines. Therefore, if we define: 



N = N, + Nk + N, 



N 



we obtain the quantity: 



iV<^ 1 + 



■ cos [(ka 



(kfc - kc) • Ti + $ 



and the probability is proportional to 

p 



d$ r'^ rf* -TT f 3 



cos 



(kfa - kj • r, + $ 



(kc - ka) • r, - $ - $ 



cos (kc - ka) • r, + -$ - $ 



(17) 



(18) 



(19) 
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The situation is therefore a direct generaHzation of that studied in the preceding section. Initially, the two phases $ 
and $ are completely undetermined but, as more and more position measurements are recorded, the relative phases 
become better and better determined. We can generalize to K systems in highly populated Fock states, for which 
position measurements progressively determine the value of K ~ 1 phase differences. The phases have a property of 
transitivity (23t . .24] : measuring the a—b phase difference and the 6 — c phase difference provides by sum the knowledge 
of the a — c phase difference. 

Note however that, in all our calculations, we have ignored the intrinsic evolution of the system; as mentioned above, 
we have assumed that all measurements are made in a time that is sufficiently short to justify this approximation. 
Over longer periods of time, the inherent evolution of the condensates takes place and introduces phase dynamics. Any 
term in the Hamiltonian containing occupation numbers will tend to compete with the definition of phase introduced 
by the measurements |l5l l25l |: for instance, the mean-field interactions cause a quantum spreading of the phase |26| . 
These effects are not included in the present article. 



II. INTRINSIC AND INDUCED PHASE 



Until now, we have limited our study to initial states of the system that are products of Fock states, where the 
initial occupation numbers are perfectly well defined. We will now allow these numbers to fluctuate, beginning with 
the simple case where their difference fluctuates but their sum remains constant. 



A. Constant total number of particles 



We now assume that: 



$0 > = ^ I A^a + Q : ka; TVf, - Q : kf, > 

Q=0 



(20) 



where the difference between the population numbers varies between Na — Nf, and Na — Nb + 2Q.; we will assume 
that Q ^ Na, Nb. The calculation remains similar to that of SjJ the main difference being that a double sum: 



Q,Q' 



(21) 



is added in the calculation of the probability. In each term of this sum, instead of coming back to the same value of 
the difference Na — Nb under the effect of the sequence of a and af operators, a change Q — Q has to be accumulated 
in order to obtain a non vanishing contribution. Equation (jllf) is therefore now replaced by: 



Q,Q' {Eq^=Q-Q'} 



(22) 



Q,Q' 



2^ 



e»(Q-Q')* + e^*F+i(r,) + e-^*F_i(rO} 



i=l 



This naturally leads to the introduction of the real positive function: 



(23) 



G($) = c 



E 



XQ e 



(24) 



where c is a normalization factor ensuring that the integral of G{^) between and 2tt is 1. We can now write the 
probability as: 



r ^G(<i>)n {1 + X COS [(k, - k,) • r, + $]} 

^0 



(25) 
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which is straightforward to interpret: instead of being completely unknown between and 27r, the initial phase $ 
before the first measurements now has a distribution function G($). This function is given by a Fourier transform of 
the probability amplitudes associated with the initial variations of the difference Na — Ni, contained in l|20(l . When 
successive position measurements begin, they all contribute to a further determination of the phase <&, as in the simple 
case already discussed inQ] The information on the phase is therefore a combination of the information contained in 
the initial state and of that created by the successive position measurements. 

A classical example is a "coherent state" obtained when the coefficients xq are given by: 



XQ = e 



\a\^/2 ^ 



(26) 



where a is a parameter, which can be written: 



a = \a\ e 
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(27) 



The initial $ distribution is then given by: 



E 



„iQ(<i>-e) 



(28) 



which is a narrow distribution around $ = when \a\ ^ 1 ; in this case, the initial value of $ is already very well 
defined, so that not much additional information can be gained in successive interference measurements. 



B. General state 



In the most general case, both the sum and difference of Na and Nb vary in the initial state: 

I $0 > = E E^^» ■ N,\Na: ka; Nb-.kb > 

N Q 



where, under the sum, Na and Nh are given by: 

N 



N 



(29) 



(30) 



In the calculation of the probability, the variables N and Q play a different role, because the sequence of a and a] 
operators acts on Q but not on N . As a consequence, a double sum over Q and Q still appears, as in the preceding 
section, but only a single sum over A^. For each value of A^, the calculation is the same as above, and leads to the 
definition of a A'^ dependent distribution function for $: 



Q 



(31) 



Then the sum over A^ can be made in a second step; one can easily see that equation l|25|l remains valid provided one 
introduces the following definition of G(<I>): 



(32) 



N 



We therefore find a situation where each value of A^ contributes independently to the distribution of the phase $. 
Different situations are possible. If for instance all functions Gn{^) are peaked around the same value of the phase, 
we obtain a state where the initial phase is well defined, as for the coherent state (|26|l . But if, on the contrary, they 
are either non-peaked functions, or functions peaked around different values, the initial phase is uncertain; it only 
becomes better and better known as more and more position measurements are performed, as in 

This shows 

that the spontaneous appearance of relative phase that we have discussed is not restricted to Fock states; it is actually 
a general property of all states involving high population numbers, provided the number of measurements remains 
much smaller than these populations. 
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III. SPIN STATES; QUANTUM NON-LOCALITY 

We now apply our calculation to particles with two internal states, which we note a and (3. We consider that these 
states are spin 1/2 states; the notion of pseudo spin allows us to do so without loss of generality. 

A. Calculation of the probability 

The initial state of the system is assumed to be: 

I $0 > = \Na: ka,a ; Nt : kf,,/? > (33) 

We note ^/j(r), with ^ — a,/3, the field operators associated with internal states a,/3. The r dependent local density 
is then: 

n(r)= *]j(r)^'^(r) + *l(r)«'„(r) (34) 
and the three components of the local spin density are: 

a,{r)= v|/t(r)v|/^(r)-vl/t^(r)vl/„(r) 

<Jx{r)= ^'j3(r)^'„(r) + «'l(r)*;3(r) ^ (35) 



f^yW- i *]5(r)«'a(r)-«'t,(r)*;3(r) 



The spin component in the direction of plane xOy making an angle 6 with Ox is: 

ae{r) - e^' ^l^r)^ ^{v) + e-''^Ur)M' p{r) (36) 
This operator, in spin space, has two eigenvalues 77 = ±1; the corresponding projectors are: 

P,(r,0)= i[n(r)+,/ae(r)] (37) 

We now consider a measurement of one spin at point ri in direction 0i, another at point r2 in direction 02, etc.; 
this corresponds to the sequence of operators: 

o-ei(ri) ; ag^{r2) ; (T93 (rg)... crgp (rp) (38) 

As in m we assume assume that all r's are different so that all these operators commute; in addition, as already 
mentioned, in all this article we assume that the sequence of measurements is sufficiently brief to ignore any intrinsic 
evolution of the condensates. Using again the expansion of the field operators on the annihilation operators: 

^^(r) ^ ^e*-ak;^ (39) 

k 

we obtain the probability of a sequence of results: 

Tyi = ±1 ; ?72 = ±1 ; ... VP = ±1 (40) 
in the form of an average value in state | $0 > ^ 



< $0 I P^i(ri,6ii) X P^,(r2,6i2) x x P^^{rp,9p) \ $0 > 



e * a, , a. 



k';a"'k.;/3 



$0 > 



(41) 



The rest of the calculation is now very similar to that of ijl Al In each term contained in the product of the second 
line of H41|l. we may assume that the a and a' are ordered in the "normal" order (with all the a's operators to the right, 
all a^'s to the left). This is because, in the first line of H41(l . when the projectors are replaced by their expressions in 
function of the field operators (and hermitian conjugate), one can move all ^f's to the right, all \l/^'s to the left (since 
operators at different points of space commute). Now we see that, each time a k^ or a kj appears associated with 
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the spin index a (or (3), it is necessarily equal to (or k^) to give an non-zero contribution; the operators a and 
introduce numbers y^iVajri^ as before, the only condition being that the creation and annihilation operators should 
be balanced in each sequence. If we assume as above that the number of measurements P is much smaller than N, 
the second line becomes proportional to: 



/'^ ? n [l + I n^ e'[('^"-'''')-+'^'l e-* + C.C.)] (42) 
•^0 i=l 



where the contrast ratio x is defined by equation H1U|) . The probability of the sequence of results rji, r]2, ..rip obtained 
at points ri, r2 , ..rp is then: 



/■^'T ^(J) 

/ {1 + a; 77i cos [(ka - kfc) • r, + 0i - $]} 

Jo 2^ t\ 



(43) 



which is the general result. 

From now on in this section, we assume that Na = Ni,, so that the contrast ratio x is equal to 1. For a sequence of 
+1 results, iPHjl reduces to: 



TT 2 

t;— cos^ 



(k, - kb) • n+f?, - $ 



(44) 



A simpler case occurs when ka = k;,. The probability of a sequence of P+ results equal to +1 and P_ results equal 
to —1 then becomes: 



2^ 



llcos (^— )_[|sm (^— ) 



(45) 



or, if all angles 6 are equal: 





$■ 


2P+ 






cos- 






lo 2^ 






sm-^ 



1 2P- 



(46) 



(Walhs integral). 

If we compare with the situation studied in ^ we now have an additional element: the adjustable parameter 9i 
(direction of the measured spin component in the xOy plane) . This combines with the dependent interference effect 
to determine the probability of each sequence of results, as shown by equation 1031) • When ka — kb, these spatial 
interference effects disappear, and the probability of the succession of results takes the simpler form (|45|l . where the 
effect of the t^'s and 6''s is isolated. For the first measurement, the two results rj — ±1 are equally probable but, as 
soon as one of these results is obtained, the distribution of $ is changed. The basic process behind the progressive 
determination of $ is the same as in ^ after a given sequence of measurements with known results, the new 
distribution of probability for $ is given (within a normalization constant) by the initial probability of this sequence, 
i.e. by replacing in (|43|l the measured ry^'s by the corresponding results. Each result rji = +1 brings in a factor 
cos^(6'j — <I>)/2 that tends to localize the $ distribution around ^ — 9i and makes it vanish in the opposite direction 
^ = 9i + tt; each result rji = —1 brings in a factor sin^(6'i — $)/2 that does exactly the opposite. A combination of 
two successive opposite results rj = +1,-1 brings in a factor ~ sm^{9i — $) that cancels the $ distribution in the 
two directions $ = 0^ and ^ = 9i + tt and tends to localize it around ^ — 9izL 7r/2. Similarly, sequences of rjiS with 
different numbers of results +1,-1 will introduce a $ maximum in some intermediate direction (see for instance the 
discussion of §III-C of ref. 

This illustrates how a sequence of successive measurements can progressively make the $ distribution more and 
more peaked around some value $ = 'I'max, and brings the system closer to a situation where all the spins are polarized 
in the same transverse direction 9 — ^mstx-', interesting numerical simulations of the phenomenon can be found in 
reference 27j. The individual discrete measurements cooperate to provide information on a quasi-continuous quantity, 
the direction of the transverse orientation. Moreover, the experimenter can adjust the choice of the measurement angles 
9i in function of the preceding results ±1 obtained. It is not necessary to orient the measurement apparatus in the 
direction of the component that is to be measured, as would be the case for a single spin; a possible strategy to 
measure the position of the maximum $niax is to perform measurements in the perpendicular direction and check 
that the proportion of 77 ± 1 results is 1/2 (in this way, first order 9 variations can be measured). 
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B. Quantum amplification 

An interesting property of the spin measurement sequence is that it can give information on the transverse spin 
orientation of a very large number 2N of particles even when P, the number of actually measured particles, is much 
lower. In principle, one can get information on the spin orientation of a macroscopic sample by just measuring the spin 
orientation of a microscopic sub-sample made of 100 particles for instance, which involves a huge amplification factor. 
Such an amplification is often evoked in the theoretical discussion of the measurement process in quantum mechanics, 
where a microscopic system (the measured system) triggers an instability in the measurement apparatus, leading to 
macroscopically different states of the pointer (here the large number of non- measured particles). The phenomenon 
is clearly related to the effect discussed by Siggia and Ruckenstein '28!| in the context of the spontaneous appearance 
of a transverse spin orientation in a double condensates of spin polarized atomic hydrogen. 

A related discussion is given by Leggett and Sols in ref. [l^ (see paragraph before last of this reference), who 
consider two superconducting systems that are initially described by an incoherent mixture of number states, and 
between which a Josephson current can flow. They then ask the question "does the act of looking to see whether 
a current flows itself force the system into an eigenstate of current and hence of relative phase?" . They assume the 
presence of a small compass needle measuring the magnetic field produced by the Josephson current. Both the needle 
and the current are macroscopic, but one can assume that the current is arbitrarily large while the needle is tiny. 
Under these conditions, it seems "bizarre in the extreme" to assume that, by some mysterious amplification effect, 
it is the small classical object that will force the large one to take a definite value. The authors conclude that, while 
standard mechanics would answer "yes" to their question, "common sense rebels against this conclusion, ... and 
we believe that in this case common sense is right" - in other words they conclude that the phase existed before 
the measurement (additional variable). In our case, the discussion is different since no external classical pointer (the 
compass needle) is introduced. In a sense, the point is even stronger since the small object forcing the large classical 
object into a definite value may even be microscopic. But the conclusion remains the same: if one rejects the idea of a 
small object creating the value of a macroscopic variable (the phase), then one is led naturally to accept the existence 
of hidden variables in quantum theory. We note in passing that, as remarked by Bell pgllsof . even if it is traditional, 
it is somewhat clumsy to call "hidden" an additional variable such as $: it is precisely the variable that is actually 
seen in the experiment, while the variables of quantum mechanics (wave functions) tend to remain invisible. 



C. Non-local quantum effects 



For the discussion of this subsection, it is convenient to assume that the orbital states of the particles are not 
necessarily plane waves. We therefore replace (|33ll by: 

|$o> = \Na:^a,a;Nb:^b,l3> (47) 

where the orbital states correspond to the wave functions: 

< r I .^a >= '^a(r) , < r I ^f, ifb{v) (48) 

with relative phase given by: 

arg{^a(r)/^6(r)} =^(r) (49) 

The calculation remains similar to that given above. In H39|) . the field operator was expanded on the annihilation 
operators at corresponding to plane waves; here we expand (r) on annihilation operators corresponding to a base 
where the particle is in internal state a and in orbital states with wave functions 'y3a(r), <P2(r), </53(r), etc.: 

oo 

- ^'/'^(l-) ag; M (50) 
q=l 

with 'Pi(r) ~ ¥'a(r), and where 932(1"), V?3(r), etc. complete V'a(r) in order to make an orthonormal basis in orbital 
space. We do something similar for ^'/3(r), now in a basis with orbital states that contain (y5b(r) as the first vector. In 
practice, the only difference is that the exponentials g^A'Mo.-'i^b)-ri\ j^q.^^ replaced by the product 93^(r)(pa(r)- The 
result for the probability is then: 

/•2^ (f$ r 1 

- —X{\^^a{v^)\' + Wb{v,)\' + 2xl^,\^aiy,)^b{v^)\ COS [^(r,) - 0. + $]J (51) 
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Even if Na = (so that x = 1) the contrast ratio of the interference term, which contains the 6i dependence of the 
result, is not necessarily 1; it is maximum at points r where they have the same modulus, zero if one of them vanishes. 
Clearly, in order to determine the relative phase $, transverse spin measurements should be made in the regions of 
good overlap between ipa (r) and ipt (r) , 

We can now assume that fa(j) and (/'{)(r) are equal and that both these functions arc non-zero only in two distant 
regions of space, D and D ; for simplicity we will assume that they are constant within these two domains; equation 
(|51|l then simplifies and the orbital wave functions disappear, as in (|45|) . All the measurements are made in region D, 
which is supposed to be relatively small so that it contains a small average number of particles, 10 or 100 for instance. 
Region D can be much larger and contain a macroscopic number of particles. Both regions can be very distant 
from each other, so that no signal transmission between them is possible during the duration of the experiment. The 
curious prediction of quantum mechanics is that a measurement of a microscopic transverse spin orientation in region 
D will immediately induce the appearance of a parallel macroscopic spin orientation in region D . 

One can make the situation even more striking by assuming that the wave functions (/?a,b(r) are non-zero in three 
regions of space: a small measurement region D as before, a large region D at a, small distance from D, and finally 
another large region D very far away (in a different galaxy for instance). The transverse spin orientation in D then 
plays the role of a local pointer, which can be seen as a part of the measurement apparatus which measures the 
direction of the transverse spin orientation in D. Remote region D contains another pointer, which under the effect 
of a measurement in D immediately takes a direction parallel to that of the pointer in D ; the phenomenon could then 
be called teleportation of the direction of pointers, a clear illustration of quantum non-locality. Of course, we should 
also make the usual proviso: the direction obtained in the measurement can not be controlled, but is completely 
random: there is no way to use the teleportation to transmit information UM, a necessary condition for preserving 
consistency with relativity. 

D. The EPR argument for condensates 

The perfect correlation between the direction of pointers at arbitrary distances is of course reminiscent of the 
quantum correlations discussed in the context of the famous Einstein-Podolsky- Rosen (EPR) argument; the similarity 
is even more striking in the EPRB version proposed by Bohm (B), which involves the measurement of spin directions. 
For a review, see for instance ref. and references therein; a detailed historical perspective is given in 33] . EPRB 
consider two correlated particles which undergo spin measurements in two remote regions of space and assume that 
quantum mechanics gives correct predictions concerning the probabilities of the various measurements of the spins 
along different directions. The EPR reasoning [s^ states that "..if, without in any way disturbing a system we can 
predict with certainty the value of a physical quantity, then there exists an element of physical reality corresponding to 
this quantity" (may the most quoted sentence of all physics!). Since EPR attribute independent elements of physical 
reality to the content of remote regions of space, this is often called an assumption of "local realism" . From this they 
show that some elements of physical reality are not contained in quantum mechanics, in other words that this theory 
is incomplete. Bohr rejected this conclusion 35] because he considered that the notion of element of reality used 
by EPR contained essential ambiguities. In his view, "the procedure of measurement has an essential influence on 
which the very definition of the physical quantities in question rests" . One should then only ascribe physical reality 
to the whole extended system made of the microscopic particles and the macroscopic measurement apparatus, not 
to microscopic subsystems only. As a consequence, one can not always attribute distinct physical properties to the 
content of in different regions of space, as EPR do. In other words, Bohr sees the process of measurement on an 
extended microscopic system as a fundamentally random process that is not necessarily localized in space js^ . 

Let us now apply the EPR reasoning to the evolution of the macroscopic pointers (transverse spin orientation) in 
D' and D , when a series of spin measurements is performed in D. We focus for a moment on D and consider 
the elements of physical reality associated with this region of space. As EPR, we consider that they can not vary 
suddenly under the effect of events taking place at arbitrarily large distances in D and D . Now compare the elements 
of reality associated to region D , just before, and just after the series of spin measurements performed in region 
D. After the measurement the region contains a macroscopic spin orientation in a given transverse direction, which 
are necessarily associated with some elements of reality. But these elements cannot appear as a result of a random 
process taking place in Z?, with no possible causal link since the distance is much too large. Therefore, even before the 
measurement, the element(s) of reality associated with the spin orientation already existed. But, on the other hand, 
standard quantum mechanics does not include this element before the measurement, only after; it therefore misses at 
least one element of reality, it is incomplete. 

One can make the reasoning more explicit by considering two realizations of the experiment, one in which the 
transverse spin orientation is found at the end of the experiment in a direction $ = $1 (with some small uncertainty 
A<I> depending on the number or measurements), another in which this spin orientation is found in a different direction 
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<I> = $2, after exactly the same sequence of measurements. Standard quantum mechanics describes the two final 
situations by different vectors in the space of states of the system, accounting for the existence of different physical 
quantities in region D (different elements of reality). On the other hand, it considers that the initial states were 
exactly the same; the differences appear only under the effect of the measurements. But this is in opposition with the 
EPR notion of local reality, which implies that the differences must have already existed before the measurements. 
In other words, the EPR reasoning shows that the two realizations of the experiment actually started from different 
initial conditions. Since quantum mechanics ignores this difference, it is incomplete. 

What are the consequences of this transposition of the EPR reasoning to a different situation than the two par- 
ticles of the original article? First, the EPR sentence quoted in the first paragraph of this subsection becomes even 
more true. Here the pointers in D and D are not invisible microscopic objects, but macroscopic and permanent 
angular momenta, similar to the spin magnetization of ferromagnets; they are classical physical entities that can be 
manipulated directly, so that it seems difficult to deprive them of properties that are independent of the measurement 
apparatuses. Therefore Bohr's denial of independent physical reality becomes more difficult to accept (of course, no 
one knows if Bohr would have given again it in this context!). Second, even if one follows Bohr and takes for granted 
that physical properties are reserved to systems including macroscopic measurement apparatuses, another difference 
arises here: one single experimental setup can lead to many spin orientations. With a single spin, one can measure 
only the spin component along the direction defined by the apparatus; with many spins in region D, it is the measure- 
ment process which determines the axis along which the spins become oriented (the transverse orientation can take 
in general a direction which makes any angle with the direction of measurement). In order to show that incompatible 
observables can simultaneously be defined (as opposed to standard quantum mechanics), the usual EPR argument 
involves different (and incompatible) experimental setups. Here, no measurement at all is performed in D and D ; 
moreover, the sequence of measurements performed in D may be always exactly the same in successive realizations of 
the experiment, even if they lead to different directions of the transverse orientation Clearly, Bohr's argument 
ascribing different physical realities to incompatible measurements does not apply anymore when only one kind of 
measurement takes place. Another difference is that, in the usual discussion of EPR correlations, the entanglement 
of the two spin particles plays an essential role; here, the initial state vector is a simple product involving the two 
internal spin states. One can finally remark that, in usual EPRB situations, one has to invoke a "no quantum cloning" 
theorem to explain why superluminal signalling is impossible. Here, in all regions D, D and D , we have a large 
number of copies of the individual systems to measure precisely the transverse direction of polarization. 

The EPR argument therefore certainly looks different, and probably even stronger, in the context of Bose-Einstein 
spin condensates. Within standard quantum mechanics, it seems really difficult not to speak of action at a distance 
to describe the instantaneous appearance of the macroscopic transverse orientation of the pointer in region D - but, 
again, we must remember that its final direction is totally uncontrolled and cannot be used as a signal. In other 
words, what we we have is an "uncontrolled action at a distance". On the other hand, we do not have the equivalent 
of a Bell theorem in this case: the introduction of an additional variable $ to quantum mechanics easily allows one 
the reproduce the prediction of quantum mechanics in a purely local model. Amusingly, the situation is exactly the 
opposite of that usually discussed in the context of EPR-Bell experiments on pairs of spins: here non-locality appears 
clearly in standard quantum mechanics, but disappears in theories with additional variables. 

E. Conservation of angular momentum 

An intriguing question relates to the momentum absorbed by the measurement apparatus during the interaction 
with the measured system. For instance, in a usual Stern-Gerlach experiment where a single spin, initially polarized in 
a transverse direction, is measured to end up in a longitudinal direction (with respect to the measurement apparatus) , 
it is generally assumed that the variation of its angular momentum is absorbed by the apparatus - this variation is 
microscopic and so small that it is impossible to measure in practice. Here, we have a system which starts before 
measurement from a state ()47|) with zero average value of the angular momentum 0; after a sequence of measurements 
performed in D it reaches another state where all spins point in a transverse direction given by the spontaneous value 
of which has a high angular momentum if there are many spins. From the usual conservation rule, one could then 
expect a priori that the measurement apparatus should acquire the opposite amount of angular momentum. On the 
other hand, we have seen that the measurement of a few spins in D may result in the appearance of a macroscopic 
angular momentum for the whole system in D + D . Does the measurement apparatus really absorbs all the angular 
momentum acquired by the whole system, although it interacts directly only with a subsystem in D that is very 
small? 

The usual rules of quantum mechanics are not very specific about the answer to this question. The emphasis is 
usually put on measurements of microscopic systems, which have a completely negligible effect on the macroscopic 
measurement apparatus - an exception, nevertheless, is the discussion of the recoil effect of a moving double slit ex- 
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periment in the famous Bohr-Einstein debate. Here, due to the amphfication process discussed in Wl Bl the measured 
system itself is macroscopic. If we assume that the answer to the question is "yes" , in order to keep a strict conser- 
vation rule for the angular momentum, we open the way to paradoxes: for instance, by locally dephasing in region 
D the two wave functions ipa and ipt, before the measurement in D, one would control the direction of the angular 
momentum in D induced by this measurement, and therefore the transfer of angular momentum to the apparatus 
in region D. This would allow superluminal signalling, in contradiction with relativity. The appropriate answer is 
therefore probably rather "no, the angular momentum transferred is just that of the individual spins that are mea- 
sured" ; in this perspective, the measurement operation creates in region D an angular momentum that appears from 
nothing, which is somewhat paradoxical too. Maybe this just means that the rule of angular momentum conservation 
should to not be taken too strictly in a quantum measurement processes. Another possibility, again, is to take the 
point of view of theories with additional variables, and assume that $ was initially actually perfectly determined in 
each realization of the experiment, even if it is unknown; in this perspective no angular momentum transfer has to 
take place, so that the difficulty vanishes. 

F. More than two internal levels 

The above considerations can easily be transposed to more than two internal states, three for instance. The situation 
then becomes analogous to that of 31 Bl where two relative phases $ and $ enter the calculation instead of one. One 
can assume that transverse spin observables relative to the a — /3, and to the /? — 7, transitions are measured. 
A sufficient number of measurements will determine the two relative phases. An interesting property is that the 
transverse orientation of the spin associated to the a — 7 transition will then be determined, although no direct 
measurement of this observable has been made. 



IV. CONCLUSION 

In standard quantum mechanics, physical systems in highly populated Fock states have a completely undetermined 
relative phase, but this phase can become very well defined under the effect of a few interference measurements only; 
this is true even in the number of particles detected is negligible with respect to the total populations. The results 
mimic exactly what would happen if the initial phase was fixed, but initially unknown, and progressively revealed 
by the observation of the interference pattern. So, the answer to the question mentioned in the introduction "do 
condensates that have never seen each other have a relative phase?" is: in standard quantum mechanics, in principle 
they do not, but they can acquire it under the effect of a few microscopic measurements, in a way which gives the 
impression of a well defined initial phase. As Anderson puts it T^: "any future experiment can be interpreted as 
if $ was fixed" . It is therefore tempting to assume that the phase existed from the beginning. This is the point of 
view of physicists who argue that, for superfluids, the usual postulates of quantum mechanics should be completed by 
an additional postulate introducing spontaneous symmetry breaking (of the conservation of the number of particles) ; 
this leads to non-zero average value of the field operator. This view is of course legitimate, but is actually nothing 
but another form of the "hidden variable" theory of quantum mechanics (de Broglie, Bohm), as we have discussed 
m § It also remains legitimate to stick to orthodox quantum mechanics and to consider that it is the process of 
quantum measurement itself that acts on the condensates and forces them to acquire a well defined relative phase. 

Spin adds to this problem the notion of angular momentum and offers interesting variants of the usual EPR non- 
locality experiments. One can suppose that spin orientation measurements are performed in two (or more) very remote 
regions of space. An interesting possibility is that some of the remote systems may be macroscopic; the results in 
terms of transverse direction of the local spin orientation of the sample are not discrete, as in usual EPR experiments, 
but continuous. Actually the macroscopic spin orientation plays here the role of the usual pointer often considered in 
the theory of quantum measurement, but with a surprising property: the pointer can be put at an arbitrary distance 
from the measured system. In one region of space, a series of measurements determines a transverse orientation, in 
another remote region a classical macroscopic pointer immediately takes the corresponding direction and stays there, 
accessible to future measurements! It seems difficult to deny that standard quantum mechanics involves action at a 
distance in this case. This action is instantaneous, but not controllable, so that it can not be used to send superluminal 
signals. There has been a long debate [33, IH, H^, Eol to decide whether or not non-locality is an inherent property 
of quantum mechanics, or if it appears only if additional ingredients are added to the theory (hidden variables, local 
reality, countrafactuality [4H . etc.). Our example can be discussed just in terms of standard quantum mechanics, 
without adding such additional elements and/or incompatible experimental setups, which would allow to invoke 
different physical realities in a Bohr type argument. It therefore clearly speaks in favor of accepting non- locality as 
an intrinsic property of standard quantum mechanics. 
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